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ABSTRACT 

Context. Observation of coronal EUV spectral lines sensitive to different temperatures offers an opportunity to evaluate 
the thermal structure and flows in flaring atmospheres. This, in turn, can be used to estimate the partitioning between 
the thermal and kinetic energies released in flares. 

Aims. Our aim is to forward-model large-scale (50-10000 km) velocity distributions in order to interpret non-thermal 
broadening of different spectral EUV lines observed in flares. The developed models allow us to understand the origin 
of the observed spectral line shifts and broadening, and link these features to particular physical phenomena in flaring 
atmospheres. 

Methods. We use ideal MHD to derive unstable twisted magnetic fluxtube configurations in a gravitationally-stratified 
atmosphere. The evolution of these twisted fluxtubes is followed using resistive MHD, with anomalous resistivity 
depending on the local density and temperature. The model also takes into account the thermal conduction and 
radiative losses in the continuum. The model allows us to evaluate average velocities and velocity dispersions, which 
would be interpreted as “non-thermal” velocities in observations, at different temperatures for different parts of the 
models. 

Results. Our models show qualitative and quantitative agreement with observations. Thus, the line-of-sight (LOS) 
velocity dispersions demonstrate substantial correlation with the temperature, increasing from about 20-30 km/s around 
1 MK to about 200-400 km/s near 10-20 MK. The average LOS velocities also correlate with velocity dispersions, 
although they demonstrate a very strong scattering, compared to observations. We also note that near foot-points the 
velocity dispersions across the magnetic field are systematically lower than that along the field. We conclude, that 
the correlation between the flow velocities, velocity dispersions and temperatures are likely to indicate that the same 
heating mechanism is responsible for heating the plasma, its turbulisation and expansion/evaporation. 
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1. Introduction 

Plasma motions in the flaring solar corona are closely re¬ 
lated to heating processes and can be used to study the 
spatial distribution of energy release and dynamics of en¬ 
ergy transfer. New instruments onboard recent solar space 
missions, such as Hinode/EIS and SDO/AIA, offer an op¬ 
portunity to study characteristics of very hot (~1-10 MK 
and hotter) pl asma in flaring atmosphere with high spa - 
tial resolution ( Culhane et al.[ 2007 [Lement et al. 20121. 
Observations show that solar flares have a complex plasma 


flow structure at different spatial scales (e.g. Watanabe et 


al. 

2010 

Del Zanna et al. 

2011 


containing both red-shifted and blue-shifted components 


(e.g. Young et al.[ |2013[ ), indicating that plasma motions 
are inhomogeneous in spatially unresolved volumes. 

Analysis of the Doppler shift structure in coronal lines 
shows that colder plasma (T « 0.1 — 1 MK) normally moves 
downwards with velocities of about 20-50 km s“^, while the 
hotter plasma (T « 1—10 MK) shows upflow; this upflowing 


plasma is very inhomogeneous, showing vel ocities from only 
few km s“^ up to 400 — 500 km s“^ (see Doschek et al. 


aTTpOlO 

^o not& 
denser 


0 


Dosche 
that t. 


2007[ Del Zannaj 2008[ Doschek et al.[|2008 Watanabe et 


K et al^2013 Young et al. 2013^. It is 

le downfiow plasma can be substantially 
n ~ 3 X 10^® m“^) than the upflow plasma (n ~ 
5 X 10^“^ — 10^® which is, probably, why the downflow 

is easier to observe than the upflow (e.g. Del Zanna| 2008 
Doschek et al.l|2008|. 
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Most coronal lines in flares also demonstrate non- 
ther mal broadening due t o additional velocity dispersion 
(e.g. Doschek et al. 2007| . The non-thermal velocity ap¬ 
pears to~correIatewitE~the bul k plasma velocity: for in¬ 
stance, in two flares analysed by [Doschek et al. (20081 us¬ 
ing Hinode data, non-thermal velocity dispersion is higher 
at locations with higher upflow speed: the velocity disper¬ 
sion is about 25-30 km s“^ when the upflow speed is nearly 
zero, and is about 50-60 km s“^ when the upflow veloc¬ 
ity is 20 km s“^. The non-t hermal velocities, w hich appear 
early in solar flares (see e.g. Harra et al. 2013), have been 
found to correlate also with plasma temperature. The non- 
thermal velocities have lowest values near 1 MK: from 0 to 
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Gordovskyy et al.: Plasma motions in twisted loops 



Fig. 1. Evolution of Model LI. Left column: Magnetic field lines; red and blue field lines are tied to opposite footpoints 
of the twisted loop, yellow lines denote ambient field. Middle column: Temperature distribution; yellow patchy surface 
corresponds to 2 MK, orange and blue surfaces correspond to 3 and 4 MK, respectively. Right column: Surfaces of 
j = jcrit^ i.e. showing the locations with fast magnetic energy release. Times (after onset of the instability) are shown 
for each row on the left in units of characteristic Alfven times tA, which are 0.44 s for this model. 


20-40 km s“^. At lower temperatures (0.1 —1-MK) these ve¬ 
locities are higher, about 40-60 km s“^. At the same time, 
at high temperatures (10-20 MK) they can be as much as 
100-130 km s“^. The non-thermal velocities at high tem¬ 
perature vary during the flare, they peak very close to the 
hard X-ray maximum at about 100 km s“^ (some o bservers 


note much higher values, up to 380 km s see e.g. Cirtain 


et al. (20131), and then duri ng the decay phase of the flare 
reduce to about 50 km s“^ (Susino et al., 20131. 


The bulk plasma flows are normally interpreted as a 
vertical advection of plasmas in stratified atmosphere due 
to heating in the chromosphere, transition region or lower 
corona. As far as the non-thermal line broadening is con¬ 
cerned, there are several viable interpretations. Firstly, this 
can be due to the turbulence, which can occur in a very 


hot plasma in post-reconnection loops or directly at the 
primary energy release location. The latter is compatible 
with the idea that the reconnection should occur in a very 
turbulent medium, and this t urbulence would facilitate the 
magnetic diffusion (see e.g 


Lazarian & Vishniak, 1999 


Browning & Lazarian 2013). It is also expected that the 


plasma in flaring coronal loops will be strongly turbulent 
due to relatively fast heating (by energetic particles, waves, 
shocks, conduction etc). Secondly, coronal lines can also 
be broadened due to unresolved inhomogeneity of regular 
flows (either in the plane perpendicular to LOS, or along 
LOS). This is possible, particularly, in the lower corona and 
the chromosphere, which can have small-scale structure (as 
small a s ^100 km, see e.g.|Antolin fc Rouppe van der Voort| 
(2012); Gordovskyy & Lozitsky (2014[)). The velocity field 
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can also be highly non-uniform in a relatively small volume 


Vrsnak et al. 

2009 

Ko et al. 

2010 


(e-g- _ ^ _ 

The issue of plasma turbulence in solar flares is inti¬ 
mately connected to the problem of energy partitioning. 
There are several estimations of the energy partitioning in 
flares, but they appear to be very different, changing from 


flare to flare (see e.g. Emslie et al. 2004 Fleishman et al. 


20151. The problem is that all observational estimations ap¬ 


pea r to depend strongly on the instrument properties (see 
e.g. Benz 20081. Furthermore, there is a substantial un¬ 


certainty about the total energy of non-thermal particles: 
it depends on the low-energy part of the spectra, which, in 
turn, depends on the so-ca lled lower energy cu t-off, which is 


still under discussion (e.g. Kontar et al. 2008). Diagnostics 
of thermal plasma can also be difficult, particularly of low- 
temperature plasma, which is less visible in FUV continuum 
and contributes less to coronal spectral lines. Indeed, an ad¬ 
equate evaluation of the turbulent velocities is necessary in 
order to determine what part of the energy is released as 
either kinetic, thermal, magnetic or potential energy. The 
first two types depend on the plasma flows, turbulence, tem¬ 
perature and density distributions in the plasma, and ener¬ 
getic (i.e. non-Maxwellian) particle spectra. In this study, 
we derive the parameters of thermal plasma from our MHD 
models. 

Most studies interpret these observations based on 
the so-called standard model of solar flares. At the same 


unstable twisted coronal loops (Browning et al. 

2008 

Hood 

et al.l 

2009 

IGordovskyy &; Browning 

|2012||Baretord et al.l 

2013 

Gore 

ovskyy et al.[ 

2013 

2014 

, which are more rel- 


evant to confined hares within single loops, offer a good 
opportunity to investigate plasma motion in flaring loops 
at different spatial scales, from ~100 km to ^10 Mm. 
One advantage of this scenario is that it offers both con¬ 
figurations with rather localised heating (at the loop-top 
and/or near foot-points), as well as configurations with 
heating sources quite uniformly distributed along flaring 
loops. Characteristics of plasma motion are defined pre¬ 
dominantly by the distribution of plasma heating sources 
(due to Ohmic dissipation, high-energy particle thermal- 
isation etc) and their temporal evolution, and by thermal 
conduction and radiative energy losses. Hence, even though 
the considered models are not a universal flare scenario, and 
don’t account for energy transfer by energetic particles, the 
results obtained with these models can be extrapolated to 
a wide range of flare configurations and scenarios. 

In the present paper we study turbulent plasma mo¬ 
tions, resulting in non-thermal broadening of spectral lines. 
Our models make it possible to investigate the velocity field, 
with spatial scale of ^100 km (which is the grid-step in 
z-direction), which is not resolved in majority of observa¬ 
tions. We derive characteristics, which are direct proxies 
for measured parameters, and discuss them in context of 
observational data. 


2. Main features of the considered flare model 

We consider four different models of kink-unstable coronal 
loops with parameters typical for a flaring atmosphere: 

— Model SI has a loop lengths of about 20 Mm with 
the apex at about lOMm. Its footpoint magnetic field 



Fig. 2. Locations of sampling boxes in the domain: green 
box denotes the “footpoint region”, red box denotes the 
“coronal region”, and the blue box denotes the “loop-top 
region”. For models SI and S2 the units are 1 Mm, while 
for larger models LI and L2 the units are 4 Mm. 


is 95 G. Field geometry of this loop is very simi 


lar to those described by (Gordovskyy et al. 2014 


Pinto et al. 2016). The model box has dimensions 
X = [—10;-1-10J Mm, y = [—10;-1-10] Mnn and z = 
[0; -1-20] Mm. The axes are shown in Figure!^ 

Model S2 is considered in the same numerical box 
and has the same parameters as model SI, apart from 
stronger footpoint magnetic field - it is 200 G. 

Model LI has similar configuration to smaller loops {SI 
and S2) but is proportionally larger and has stronger 
field. Its length is about 80 Mm and footpoint strength 
is 700 G. The numerical box is x = [—40; -1-40] Mm, 
y = [—40; -1-40] Mm, and x = [0; -1-80] Mm. 

Model L2 is considered in the same numerical box 
and has the same parameters as model LI, apart from 
stronger footpoint magnetic field - it is 1500 G. 


It is assumed that in all the four models the lower bound¬ 
ary corresponds to the same level in the solar atmosphere 
and, hence, the lower boundary has the same density in 
all models, 10^° m“^ (approximately corresponding to the 
upper chromosphere). 

The loops are all embedded into gravitationally strati¬ 
fied atmosphere. The initial density in the atmosphere is 


Pt=o{z) = Pi exp 


Zscl ) 


+ p 2 exp 


Zsc2 ) 


( 1 ) 


where the density constants are pi = 3.34 x 10“® kg m“^ 
and p 2 = 2x 10“^^ kg m“^, and the density length scales are 
Zscl = 0.25 Mm and Zsc 2 = 50 Mm. The shift Zsh=^-5 Mm, 
putting the lower boundary of the computational domain 
approximately at the height of the lower chromosphere. 
The initial density at the upper boundary is 10^® m“^ 
in the models SI and S2, and 3 x 10^^ m“^ in LI and 
L2. The gravity acceleration is assumed to be constant 
g = 275 m s“^, so that the initial pressure can be easily 
obtained from density, as dpt=o{z)/dz = —pt=Q{z)g. As the 
density constants pi and p 2 , and the length scales Zgci and 
Zsc 2 are very different, there are two regions with nearly 
constant temperatures: around 10^ K below 2 Mm, and 
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about 0.9 MK about 2 Mm (see [Gordovskyy et al. 20141. 
The magnetic field is strongly convergent, i.e. it is about 
10 time weaker near the loop apex, compared to the foot- 
points. This means that the fluxtube cross-section radius 
near the apex is around 3.2 times larger than that near 
foot-points. The modelling has been performed using the 
LARE3D code by Arber et al. (20011. The simulations are 
based on the resistive single-fluid MHD, also i ncorp orat- 
ing Braginskii thermal c onduction ( Bragmslm 19651 and 
radiative losses following Klimchuk et al. ( [2008 ). 

Pre-kink configurations have been derived as described 


Gordovskyy et al. (2013); [Bareford et al. (2016); Pinto 
et al. (|2016p. Although the initial atmosphere is in the hy- 


drodynamical equilibrium, it is not stationary because of 
the thermal conduction. The heat flux from the hot corona 
to colder chromosphere results in noticeable changes of the 
temperature around the transition region. Moderate heat¬ 
ing of the transition region and chromosphere and moderate 
cooling of the lower corona lead to smoother temperature 


profile near the transition region (see Pinto et al. 2016 


Bareford et al. 2016). The atmosphere settles in about 
1000 t^ and only after this the initially potential loop is 
twisted by footpoint rotation. The twisting is very slow, the 
maximum linear velocity in small models is about 10“^ vq 
( or ~ 3 km s^), which is always lower than local Alfven and 
sound velocities. Hence, during the twisting phase loops un¬ 
dergo sequences of nearly force-free states until kink insta¬ 
bility onset. As the result, the loop undergoes a sequence of 
nearly-force-free states during the twisting. Also, the twist¬ 
ing velocity is well below the observed LOS velocities and 
velocity dispersions (> 10 km s^). 

The kink insta bility occurs when the t otal twist reaches 

for details). The 
loops atmosphere 


67r — Stt (see also Bareford et al. 2016 


magnetic reconnection in twisted corona. 


her studies (see 

Hood & Priest 1979 Baty & Heyvaerts 

1996 Browning & Van der Lindenl |2003| |Hood et al.l 


Gordovskyy & 

Browningl |2012| [Gordovskyy et al. 

2014 

Bareford et al. 

2016), while the evolution of thermal anc 

non-thermal rac 

iation from such systems is considered by 

Botha et al. (2012); Gordovskyy et al. (2013 2014); Pinto et 

al. (|2016|). Essentially, a loop experiences two different types 


lines of twisted fluxtube, resulting in twist reduction, and 
reconnection of twisted held lines with the ambient held, 
resulting in radial expansion of twisted loops (Figure]^. 
The plasma temperature and, hence, thermal emission in¬ 
tensities peak towards the end of fast reconnection, when 
the twist is considerably lower. As the result, the visible 
twist in EUV/SXR is much lower than the critical twist 


required for the kink instability (Pinto et al. 2016). 


3. Unresolved plasma motions in different 
temperature ranges 

3.1. Macroscopic velocity distributions and their relation with 
non-thermal line broadening 


Here we derive expressions relating the resolved how ve¬ 
locities and unresolved how velocity distribution to widths 


and p ositions of spectral lines (see e.g. Hubeny & Mihalas 
2014[ for exact derivations). 


Assume the model contains only fully ionised hydrogen 
plasma, so that the electron and proton particle densities 


and total numbers are equal, rip = rig = n and Np = Ng = 
N respectively. Further, assume that within an unresolved 
volume has a distribution of plasma in respect of the LOS 
velocity U|| and temperature T, /(i’ll, T) = Then, the 

prohle of a spectral line emitted by this volume of plasma 
will be 


I{X) = 


j J (a - Xu^g - ^Anne) dv\\dT, (2) 


T 


where a{T) is a temperature contribution function for the 
spectral line, s(A — Xung) is the line prohle in the absence 
of macroscopic plasma motion. 

The function s is normally a Gaussian-like prohle, which 
can be approximated by a Gaussian: 

s(A - Xiing) = So exp ^ > (3) 


where the Doppler width is AAd = y ^^Xung- (In fact, 

the exact shape of the s function is not important, as far 
as it is a Gaussian-like distribution, i.e. a distribution, be¬ 
coming 0 at ±oo, with hnite dispersion. 

Assuming that the function is also Gaussian in 

respect of the LOS velocity, i.e. 


(v) 


d?N‘^ 

dv\\dT 


= exp 


(^11 - 

^vUT) J 


MT), 


( 4 ) 


the Eq. (0 can be written (using the convolution of two 
Gaussian functions) as 


/(A) = 


t(T) exp 


T 


(A - Xiine - AXlOs)^ \ ,rp 

Axi+Axi, r"- 


( 5 ) 


Here AA„t = ^^Xune represents so-called non-thermal 

contribution to line broadening, while AXlos = ^AXune^ 
and t{T) is the product of a{T),ft{T) and constants. 

i2 7^2 

Hence, if the plasma’s distribution for each small 

temperature interval in the considered volume has some 
Gaussian-like distribution in respect of U] |, then the average 
velocity, determining the line shift can be derived as follows: 

y " 



(For this average LOS velocity we will use the notation {v) 
hereafter.) The half-width of this distribution (or FWHM), 
which determines the non-thermal broadening of spectral 
line can be approximated as follows (this expression is exact 
for a Gaussian profile) 




I 


dv\ I dT 


(T'.i^||)(«|| - {v)fdV 


/ 


d‘^N^ 

dv^ I dT 


{T,m)dV 


( 7 ) 


(For this LOS velocity dispersion (or variance) we will use 
the notation Avnt hereafter.) 
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Fig. 3. Model SI. Left column: Differential emission measure versus temperature for the whole domain. Middle column: 
Average n^-weighted velocity dispersions Avnt{T) (calculated as a/Au^((T)) for the foot-point location, for x-component 
(red solid lines) and z-component (blue dot-dashed lines) of the velocity. Right column: Green dashed line is AvntiT) for 
the whole domain; red solid lines is the Au„t(T) for the coronal part of the loop; black lines are Avnt{T) for the loop-top 
region (dot-dashed line is for y-component and dotted line is for z-component of the velocity). Corresponding times are 
shown in the left panels. Locations of the sampling regions are shown in Figure 


In our simulation we use discrete forms of Equations 

and [71 


{v)t = 


AvA, = 


where 


Wti = 


o 


<5D] 


n 

EKi (^11 i - 

O 


[v)t? <5D] 


0 , 


EKi 

n 


if Tt <Ti < Tt+i 
otherwise 


In these equations, the index t corresponds to a tem¬ 
perature interval, index i corresponds to a grid points, and 


An is an elementary volume between adjacent grid points 
(constant, as the grid is uniform in each direction). Here 
Wti is, effectively, the density of plasma squared, with tem¬ 
perature Tt and velocity uiu, and, hence, corresponds to 

dv^dT fo Equations (bl?). The integration goes over 

all the grid points in the considered volume Q. 

In the above equations, the integral 

Mt = ^[wti dn] 
n 

is, in effect, the differential emission measure for the volume 
H. (In observations it is calculated as d[n‘^{T)V{T)]/dT or 
n^{dT/dl)~^ for an individual pixel.) 


In section 3.2 the velocities {v) and Avnt are calculated 
as functions of temperature using the above equations and 
compared with observational data. 
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Fig. 4. The same as in Figure]^ but for Model S2. 
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3.2. Characteristics of the veiocity fieids obtained from 
numericai modeis 

The velocity dispersion functions Avnt{T) have been calcu¬ 
lated for each model for four regions (see Figure]^: large 
coronal region (practically, the whole domain volume above 
2 Mm in models SI and S2, and above 8 Mm in models LI 
and L2), loop-top region (a cube with the size of 3 Mm 
in SI and S2, and 12 Mm in LI and L2 containing the 
loop top), the footpoints (below 4 Mm in all four models), 
as well as for the whole domain. The values of Avnt{T) 
are calculated separately using x-, y- and z-components 
of velocity. In the footpoint region, the z- and x- (or y-) 
componets of velocity correspond to directions along and 
across the magnetic field, respectively. At the same time, it 
is practically impossible to determine the dominant mag¬ 
netic field direction in the large coronal sampling regions, 
because of the twist and fast fluxtube motions. However, 
we do this for the loop-top regions. Additionally, we calcu¬ 
late the Differential Emission Measure (DEM) functions for 
each model, in order to show the distribution of the hot ra¬ 
diating plasma. The results are shown in Figures [3]|^ The 
times used here are arbitrary, measured from re-starts of 


simulations. In each model, fast energy release begins from 
the onset of kink instability, which happens approximately 
approximately at 740 s in model SI, at 270 s in model S2, 
at 310 s in LI, and at 280 s in L2. 


The DEM functions for small models (Figs. 3]|4) can be 
divided into three parts. At low temperatures ( 0.1 MK) 
there is a peak corresponding to the “cold” chromospheric 
plasma. It is followed by a nearly flat part and a peak at 
around 0.8-1.0 MK, corresponding to the coronal plasma. 
These two parts do not show substantial changes during 
loop evolution. The third part, at T > 1 MK, corresponds 
to the loop plasma heated during the reconnection. It is 
nearly flat, with slight positive slope during the reconnec¬ 
tion phase and with slightly negative slope during the cool¬ 
ing phase. This part is limited by the peak temperature, 
which is about 2.5 MK for the model SI and about 3.5- 
4 MK for S2. DEM functions corresponding to the larger 
models (Figs.[5]|^ can be divided into two parts: bulk of the 
plasma forms a peak with a strong negative slope at a few 
MK, and then show exponential profiles up to the maximum 
temperature. The maximum temperature reached during 
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Fig. 5. The same as in Figure but for large-scale Model LI. 


the reconnection phase in model LI is about 15-16 MK, 
while in the model L2 it is about 30 MK. 

Figures [3]|^ show LOS velocity dispersion as a functions 
of temperature, defined by Equation]^ The values of I^Vnt 
change substantially between different phases: they are or¬ 
der of 10^ km s“^ during the reconnection, and drop to 
^ 10 km s“^ during the cooling phase. It can be seen that 
these values normally increase with the temperature in all 
regions in all considered models. Thus, the coronal parts of 
loops in small models demonstrate an increase in Av„t from 
10-30 km s“^ to about 100-200 km s“^ when the temper¬ 
ature increases from about 0.5-1 MK to 2-4 MK. In larger 
models, this value increases from 50-100 km s“^ at T « 2- 
3 MK to 200-400 km s“^ around 20 MK. Smaller coronal 
sampling regions show a similar picture. The situation near 
the footpoints is more complicated. Overall, it is possible to 
say that in most cases the velocity dispersion is noticeably 
higher at higher temperatures. Thus, typically, in model SI 
this velocity increases from few km s“^ below 1 MK to 
about 3-10 km s“^ at >1 MK. Similarly, in model S2, it 
also increases from few km s“^ just below 1 MK to about 
10 km s“^ at ~2 MK. In large models, this velocity is even 


higher; it increases from about 10 km s“^ at 2-3 MK to 
about 50-100 km s“^ at 10-15 MK. However, near foot- 
points the Avnt(T) distribution sometimes is flat or even 
decreasing. 

Another interesting feature revealed by the non-thermal 
velocities near foot-points is the noticeable difference in dis¬ 
persions in different directions. Thus, in the footpoints re¬ 
gions of smaller loops the Avnt(T) distributions in x- and 
z-directions are similar at the beginning of reconnection. 
However, during the fastest stage of reconnection and dur¬ 
ing the cooling phase the dispersion of Vz is 3-4 higher 
higher than the dispersion of Vx- Footpoints of larger loops 
show similar picture: both during the fast reconnection and 
during cooling, the velocity dispersions in vertical direction 
are 2-3 times higher than those in horizontal directions. 
In other words, the velocity dispersion along the magnetic 
field is higher than that across the field. 

A similar picture can be seen in the small sampling re¬ 
gion near loop tops: dispersions of Vy (direction along the 
loop) are higher than dispersions of (across the loop). 
In smaller loops this is the case at three different stages, 
however, in larger loops this contrast is small during the 
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Fig. 6. The same as in Figure but for large-scale Model L2. 


onset of reconnection, becoming very substantial (up to 3 
times) towards the end of reconnection. This also might be 
explained as the effect of the magnetic field, particularly, 
towards the end of reconnection, when the azimuthal field 
is substantially reduced. There is no noticeable difference in 
the velocity dispersions in different directions for the whole 
coronal part or for the whole domain. This is quite natu¬ 
ral, taking into account that there is no substantial mean 
magnetic field in these models. 

It should be noted that the plots in the middle column 
in Figures [3]|^ corresponding to the foot-point regions sam¬ 
ple substantially large volumes (see Figure [^. Their size is 
proportional to the domain (and loop) size and, therefore, 
in models LI and L2 they sample larger volumes, including 
not only the chromosphere, but also the transition region 
and lower corona. This, most likely, explains why Avnt{T) 


functions for foot-point and coronal regions are more simi¬ 
lar in models LI and L2, compared to models SI and S2. 

It should be noted that the non-thermal velocity disper¬ 
sions near the footpoints, perhaps, are not representative of 
the whole loop: they correspond to a cooler and more struc¬ 
tured material of the chromosphere and transition region, 
where plasma /3 is higher and can be close to 1. However, 
they are still of interest, for instance, the variation of ve¬ 
locity dispersion with direction can indicate suppression or 
enhancement of the turbulence in relatively strong mag¬ 
netic field. At the same time, the coronal part represents 
the whole loop: its Avnt{T) distributions are very similar 
to those of the whole domain. 

Figure shows the time variation of the velocity dis¬ 
persion functions in three different temperature ranges, 
0.25 MK - 1 MK, 1 MK - 4 MK, and >4 MK. There is 
practically no signal in the highest temperature band in 
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Fig. 7. Variation of velocity dispersions calculated over the whole domain in three temperature bands in the four models 
after kink instabilities. Model names are shown in panels. Black solid line is for the 0.25MK < T < IMK band, red 
dashed line is for IMK < T < 4MK, and blue dot-dashed line is for T > 4MK. (There is practically no plasma with 
temperature > 4 MK in models SI and S2.) Lower limits of each time intervale approximately correspond to the onset 
of kink instability and fast energy release. 
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Fig. 8. Average n^-weighted velocity dispersions Avnt versus current densities for the four models. Each dot represents a 
cube with the size of 2.5 Lq (2.5 Mm in smaller models and 10 Mm in larger models). In order to demonstrate correlation 
with the height, colour flags indicating the height corresponding to each dot are shown in the panel S2. 


small models SI and S2. All four models show that Au„t increases with temperature most of the time. Although, in 
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Fig. 9. Distributions of current density |j| and velocity deviations 6v(x,z) around one of the energy release regions. Top 
left plot shows the global current density distribution at a; = 0 plane and the location of the sampled region. Current 
density distribution and magnetic field lines in the selected region are shown in the top right panel. Distributions of the 
p^T-weighted velocity deviations (see text for details) are shown in the bottom left (x-component) and bottom right 
(z-component) panels. (Here the units of the current density colourscale are 3.6 x 10“^ A m“^ and the units of the 
p^T-weighted velocity colourscale are 2.8 x 10^ km s“^. Distances are in 10® m. 
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Fig. 10. Average velocity (Z-component) versus velocity dispersion for model S2 (left) and L2 (right). Each point 
corresponds to a sampling cube with side 2.25 Mm in the small model (left) and 9 Mm in the large model (right). 
Here, positive velocities correspond to plasma moving up, and vice versa (contrary to the convention used in most 
observational papers). 


model LI values of Avnt in 0.25 MK - 1 MK and 1 MK - 
4 MK bands are nearly equal towards the end of reconnec¬ 
tion, while in model L2 velocity dispersions in these bands 
are nearly equal just after the onset of reconnection. 

In order to understand the relation between the primary 
energy release and velocity dispersions we calculated Avnt 
as functions of the current density j (Figure^. Similar to 
the temperature distribution of Avnt, these function were 
calculated using the original spatial resolution of our MHD 


simulations (about 0.07 Lq). Obtained Avnt{j) functions 
did not show any noticeable spatial correlation between the 
velocity dispersions and current densities. However, when 
averaged over larger volumes (cubes with about 2.5 Lq 
size), the velocity dispersions clearly show spatial corre¬ 
lation with current densities (Figure!^. Moreover, there is 
also a correlation with height. Thus, there are two types of 
hot regions in the considered models. Coronal plasma (see 
red, yellow and green dots in Figure panel S2) shows 
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weaker increase of the velocity dispersion with the cur¬ 
rent density (Avnt ~ while the hot regions in foot- 

point show much steeper increase of Avnt with currents 
{Avnt ~ j)- The effect of spatial resolution here can be 
explained by the structure of the primary energy release 
sites. Figurej^shows the (x,z) distributions of current den¬ 
sity and the velocity deviations, as well as magnetic field 
lines around a region with high current density at the ’cur¬ 
rent shield’ formed just above the top of a twisted loop. 
Here, velocity deviations are calculated as 


SVr — \l{Vr^X^z) Vr mean) 


^p^{x,z)T{x,z) 

2 ' 
Pmean^rnean 


where index r means x, y or z, and Pmean and Tmean are the 
density and temperature averaged over the sampled region, 
respectively. This region has all basic features of reconnec¬ 
tion curr ent sheet, including th e inflow and outflow regions 
(see also Bareford et al., 2016). Importantly, plasma with 
high turbulent velocities is located in the exhaust regions, 
away from the current layer, and, hence, there is no imme¬ 
diate spatial correlation between the turbulent velocities 
and current density. However, the values of Avnt averaged 
over volumes larger than these elementary current sheets 
will show spatial correlation with current densities. 

On the other hand, relatively weak correlation and a 
substantial spread of points in the Avnt versus \j \ diagrams 
means that plasma turbulence and strong current are not 
always connected. Firsly, this is because there are several 
factors, which are not directly related to the energy release, 
can result in velocity dispersion (see Sect. 3.3). Secondly, 


the magnetic energy can be converted into kinetic energy of 
plasma turbulence well away from the reconnection regions. 
This is possible, for instance, if the plasma is heated and 
turbulised by MHD shocks generated by the magnetic re- 


connection |Bareford fc Hood (2015); Bareford et al. (2016). 

We also compare the velocity dispersion (related to line 
width) with the bulk flow velocities (related to line shift); 
the measurements for our two models, S2 and L2 are shown 
in Figure They are calculated using equations |6][^ and 
integrated over temperature. Only plasma with tempera¬ 
ture higher than 1 MK is taken into account here. Each 
point corresponds to a cube with the side of about 30 grid 
points (2.25 Mm in model S2 and 9 Mm in L2). It is found 
that the average velocities are rather low when the non- 
thermal velocity dispersion is low, (v) «50 km s“^ when 
Avnt =20 km s“^. Then, it increases with Avnt, reaching 
around {v) =80-100 km s“^ for Avnt «100 km s“^, and 
about {v) =700 km s“^ for Avnt «1000 km s“^. Both neg¬ 
ative and positive average velocities are present, and the 
dependences between Avnt and {v) can be linearly approx¬ 
imated both for upflows and downflows. Although these 
dependences are clearly visible on the graphs, the spread of 
points is very big (several km s“^). 


3.3. Comparison with observations and discussion 

Our models can reproduce the two key observational fea¬ 
tures described in Section 1. Firstly, similarly to obser¬ 
vational data, the non-thermal velocity dispersion in the 
model increases with temperature, at least in its coronal 
part (which is the most of the volume for most loops). The 
exact values depend substantially on the model (or loop) 


size. At temperatures of about 1 MK, these velocity dis¬ 
persions are about 10-50 km s“^, increasing to about 100- 
200 km s“^ at temperatures of 5-10 MK and reaching up to 
500 km s“^ in the range of 20-30 MK. Observational data 
(see Sect. 1) similarly shows an increase from 10-20 km s“^ 
at ~1 MK to few hundreds of km s“^ at about 10 MK. 
Also, similar to observational data, our models show a no¬ 
ticeable peak in Avnt{T) functions at low temperatures, 
around 0.3-0.5 MK. 

Secondly, the loops considered in our models also 
demonstrate correlation of the bulk velocity with non- 
thermal velocity dispersion: the higher the LOS flow ve¬ 
locity, the higher is the Avnt value. Both, downflows and 
upflows are present (nearly equally). However, the aver¬ 
age LOS velocities derived from our models demonstrate a 
substantial spread of about 100 km s“^, while i n the obser- 
vations the (v) distribution is more compact (Doschek et 


ah] 2008). This spread is larger in the plasma with positive 


velocities (upflow). Finally, the range of average LOS veloc¬ 
ities obtained in our simulations is considerably larger (up 
to 800 km s“^) than in observations by Doschek, mentioned 
above. However, these values are still acceptable, as such 


velocities are observed sometimes in larger loop (Cirtain et 


al. 2013). 


The correlation between temperature and velocities is 
not surprising: indeed, one would expect that plasma with 
more heating would contain more kinetic energy in form 
of large-scale turbulence and regular flows. Thus, in our 
simulation highest velocities appear in and around the en¬ 
ergy release regions (close to the loop top, and in loop 
legs closer to current concentrations near footpoints). This 
also explains the correlation between the velocity disper¬ 
sion and the average flow velocity. Indeed, both values (v) 
and Avnt are high in the hot plasma regions, relatively low 
in colder regions of the flaring loop, and quite low outside 
the loop. Generally, it means that a single energy release 
mechanism is responsible for small-scale turbulent veloci¬ 
ties and large-scale plasma motion, so that the energies of 
these two types of motion are proportional to each other. 
The reason why th e s pread in velocities obtained in our 
simulations (Figure is substantially greater than in ob¬ 
servations is less clear. This may be because the large box 
we use in our measurements includes very different regions. 
On the other hand, it may be due to the fact, that in ob¬ 
servations the values of (v) and Avnt are calculated over a 
larger temperature range than in our simulations. 

Concerning the velocity dispersion, it is important to 
note that there are several different factors causing them, 
and some of these factors may be unrelated to the energy 
release. The key factor, which is in the focus of this pa¬ 
per, is the turbulence in the very hot flaring plasma. As 
shown in the Section [3.2| this turbulence occurs in the very 
hot plasma ejected from elementary reconnecting current 
sheets. Since the reconnection is the most powerful heating 
mechanism in our models, this mechanism would be dom¬ 
inant at high temperatures (1 MK and higher). Secondly, 
there are fast chaotic motions due to the loop kinking after 
the instability. Since these velocities should be highest close 
to the loop top, loop kinking would affect the temperature 
range >0.5 MK (i.e. the temperature of the quiet coro¬ 
nal plasma and the temperature of flaring plasma. There 
is also an important implication from this mechanism: no¬ 
ticeable non-thermal velocity dispersions may appear even 


before the heating, as it is found in observations (Harra et 
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al. 2013). Thirdly, the thermal flux from the corona can 
Heat dense transition regions and the chromosphere, which 
may add to the velocity dispersion at lower temperatures 
10^-10® K via non-uniform evaporation upflows. 

Substantial difference in the velocity dispersions in dif¬ 
ferent directions can be explained assuming that the mean 
magnetic field is close to the loop direction. (Here, by the 
loop direction we mean the “skeleton line” of a loop go¬ 
ing from one footpoint to another footpoint via the loop 
top.) This would mean that the velocity dispersion along 
the mean magnetic field is higher than across it. There are 
two factors, that could lead to this. Firstly, the turbulence 
in the hot flaring plasma can be anisotropic, being stronger 
along the mean field. Secondly, since the hot plasma is mov¬ 
ing predominantly along magnetic field, any cross-field in¬ 
homogeneity of the flow velocity at small spatial scales will 
result in the line-of-sight velocity dispersion, if observed 
along the magnetic field direction. This inhomogeneity is 
quite possible because of the anisotropic thermal conduc¬ 
tion: heat is quickly redistributed along magnetic field lines 
with almost no conduction in the perpendicular direction. 

Finally, there are two features in the velocity field un¬ 
related to the instability and the energy release: foot-point 
rotation and kink-oscillations of the loops before the in¬ 
stability. Thus, slow periodic variation in the total kinetic 
energ y in these loops (see e.g. Fig.4 in Gordovskyy et al. 
2014) are likely to be due to low-amplitude fundamenta 
mode (ie. N = 1) kink oscillations. Based on the kinetic 
energy variation, their amplitude should be order of O.OOIwq 
or smaller {i.e. about 1 km s“^ or smaller). Hence, these 
oscillations, in principle, could make a small contribution to 
the velocity dispersions, particularly at low temperatures. 
The maximum linear speed due to the loop rotation is about 
0.004uo, or 10 km s“^. This motion also can contribute to 
the velocity dispersion, especially at low temperature (since 
the velocity dispersion is n^-weighed and the temperatures 
in the dense chromosphere are K. However, the 

velocity dispersions higher than ~ 10 km s“^ cannot be 
attributed to the driver. 


4. Summary 

In the present work we derive characteristics of the velocity 
field in flaring loop and compare them with observations, 
focusing on the plasma velocity dispersions, which corre¬ 
spond to non-thermal broadening of coronal EUV lines. It 
is found that our models yield velocity dispersions and av¬ 
erage line-of-sight velocities, which are in qualitative and 
quantitative agreement with those derived from observa¬ 
tional data. 

Thus, velocity dispersions Avnt , average line-of-sight ve¬ 
locities (v) and maximum temperatures depend on the flare 
size and characteristic magnetic field in the loop. Thus, 
in small loops (with length of about 20 Mm) the maxi¬ 
mum temperature varies between 2 MK (loop with foot¬ 
point magnetic field about 100 G) and 4 MK (footpoint 
field 200 G), average velocities and velocity dispersions are 
about 10-100 km s“^. In large loops (with length of about 
80 Mm) the temperature can reach 16 MK (footpoint field 
700 G) and 30 MK (footpoint field 1.5 kG), while the aver¬ 
age velocities and velocity dispersions are 50-500 km s“^. 

In the chromospheric and photospheric plasma close to 
footpoints the velocity dispersion is low, about 5-20 km s“^. 


The velocity dispersion along the field is found to be lower, 
by factor of 2-4, than that across magnetic field. 

The velocity dispersion correlates with temperature. 
Thus, it is around 10 km s“^ near T « 0.5-1.0 MK, in¬ 
creases to about 100 km s“^ at ~2 MK, and reaches about 
200-500 km s-i at T « 10-20 MK. 

The velocity dispersions show spatial correlation with 
current densities, they vary approximately as Avnt 
in the corona, and as Avnt j near footpoints. 

Velocity dispersions appear to be higher along the loop 
direction, both near the footpoints and at the loop tops. 
Typically, the ratio is about 2-4. 

Finally, the velocity dispersion correlates with average 
LOS velocity. Average velocities can be both positive and 
negative (upflow and downflow). At lower Avnt values these 
two quantities are connected approximately as (v) « Avnt, 
while at high values Avnt > 100 km s“^ this relation is 
(v) « OAAvnt for upflow and (v) « O.dAvnt for downflow. 

Therefore, based on our simulations, we can conclude 
that the correlation between the observed velocity disper¬ 
sions (derived from spectral line widths) and the temper¬ 
ature, as well as the correlation between the average LOS 
velocities (derived from spectral line shifts) and the velocity 
dispersions, indicate that a single mechanism, direct plasma 
heating during magnetic reconnection, is responsible both 
for turbulisation of plasma and bulk plasma motion. 

It is important to say that a good agreement with ob¬ 
servations doesn’t mean that all flares investigated in the 
observational studies (mentioned in Sect. [H occurred in 
twisted coronal loops. However, the fact that solar flares 
occurring in this type of magnetic configurations produces 
data similar to observations implies that the ’twisted loop’ 
configuration is a viable model for a solar flare, which can 
explain observations of {v) and Avnt using simple reconnec¬ 
tion scenario with multiple extended reconnection regions. 
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